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Abstract 

A sample of 3.73 million hadronic Z decays, recorded with the OPAL detector at LEP in the 
years 1991-95, has been used to search for a narrow resonance corresponding to the decay of 
the D* /=t (2629) meson into D* ± 7r + 7r _ . The D* + mesons are reconstructed in the decay channel 
D* + — > D°7r + with D° — > K~7r + . No evidence for a narrow D* /=t (2629) resonance is found. A 
limit on the production of narrow D* /=t (2629) in hadronic Z decays is derived: 

f(Z -> 0*^(2629)) x Br(D*'+ -> D* + tt + tt-) < 3.1 x 10~ 3 (95% C.L.) 
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1 Introduction 



Quark models based on Quantum Chromodynamics (QCD) |], predict the existence of charm 
mesons whose radial wave functions are not in their ground state. A spin doublet contain- 
ing a pseudoscalar 2 1 S and a vector 2 3 Si state is expected. These states are referred to 
as the D' and D*', respectively. Their masses are predicted to be m^> = 2.579 GeV/c 2 and 
m D »/ = 2.629 GeV/ c 2 0] (see Table [1]). The uncertainty on these masses is estimated to be of 
the order of 20 MeV/ c 2 from the difference between the values obtained in computations 
and observations [[| of the masses of orbitally excited states. 

The kinematically favoured decay mode of the D* /+ is presumed to be the direct three-body 
decay into a D* + 7r + 7r~ final stateQ. This is an S-wave decay and thus generally expected to 
be broad, although a model exists 0] which estimates the partial width of this decay mode 
to be less than 1 MeV/c 2 . Decays via an intermediate orbitally excited D meson could also 
contribute to the same final state; these decay widths are estimated to be several MeV/ c 2 
by the same model. If enough phase space is available, D* /+ decays to a D* + p state with 
p — > 7r + 7r~ might also be allowed. The decays via orbitally excited states include S-wave 



(D 



*/+ 



D^°7r + , — > D* + 7r ), and D-wave transitions (e.g. D 



*/+ 



D 



*o 



7T 



D*° 

u 2 



D* 



7T 



The decays involving a p are P-wave transitions. The higher partial waves are expected to be 
suppressed due to lack of phase space. 



state 



D 

D* 
D* 
D* 
Di 
D* 
D' 
D*' 



J 



V 2 

±1/2 

°l"/2 
1 + 
-4/2 

2+ 

^3/2 
V 2 
4/2 



predicted mass 
[ GeV/c 2 ] M 



1.875 
2.009 
2.438 
2.501 
2.414 
2.459 
2.579 
2.629 



observed properties 



mass [ GeV/c z 



} 



1.865 
2.010 

(2.461) 

2.422 
2.459 



(2.637) 



width [ MeV/c z 



< 0.131 

(290) § 

18.9 
23 

(< 15) i 



Table 1: Masses of charm mesons as predicted in Ref. 0. Measurements of the meson masses 
and width are given where available || • Values in brackets are not from Ref. |§ , but represent 
recent observations from CLEO [||] and DELPHI ||. Dq, D^, Di, and DJ; are orbital excitations 
(the D* meson is also known as D^). D' and D*' are radial excitations. 



The only experimental evidence for the existence of radially excited D mesons comes from 
the DELPHI collaboration which has published an observation of a narrow resonance decaying 
to D* + 7r+7r- 0. The observed mass of 2.637 ± 0.002(stat.) ± 0.006(syst.) GeV/c 2 is close to 
the theoretical prediction. However, the interpretation of this D resonance as D*' is controver- 
sial, since the narrow width is incompatible with most theoretical models [[?]. The DELPHI 
measurement of the width of the state is limited by their detector resolution and quoted to be 
smaller than 15 MeV/c 2 at 95% C.L. 

1 Charge conjugates will be implied throughout this paper. 
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This paper presents a search for the state reported by DELPHI in hadronic Z decays recorded 
with the OPAL detector. The search covers both the region of the resonance measured by 
DELPHI as well as any other narrow resonance in the vicinity of the predicted D* /=t mass. 
Therefore the particle being searched for will be referred to as D (2629) throughout this 
paper. 

After a description of the data and Monte Carlo samples, the selection criteria for the 
D* /+ mesons are presented. Since the production mechanism of D* /+ mesons is different in 
Z — ► bb and Z — > cc events, two different selections are used to reconstruct D* /+ mesons. One 
is optimized to select D* /+ candidates in primary b-events using vertex information, while the 
second selection uses harder cuts on the momentum of the D* /+ candidate, resulting in a sample 
enhanced in cc events. In Section |j, the selection results for data and Monte Carlo events are 
presented. Systematic checks are discussed subsequently, and the calculation of limits on the 
production rate is described in Section [7|. In the concluding Section our results are compared 
with those published by DELPHI §. 



2 The OPAL Detector and Event Simulation 



A detailed description of the OPAL detector can be found elsewhere ||. The most important 
components of the detector for this analysis are the silicon microvertex detector, the tracking 
chambers and the electromagnetic calorimeter. The microvertex detector consists of two layers 
of silicon strip detectors which provide high spatial resolution near the interaction region. The 
central jet chamber is optimized for good spatial resolution in the plane perpendicular to the 
beam axisQ. The resolution along the beam direction is improved by the z information delivered 
by the silicon microvertex detector, by a vertex drift chamber between the silicon detector and 
the main tracking chambers, and by z-chambers surrounding the main tracking chamber. The 
central detector provides precise determination of momenta of charged particles by measuring 
the curvature of their trajectories in a magnetic field of 0.435 T. The solenoid is mounted 
outside the tracking chambers but inside the electromagnetic calorimeter, which consists of 
approximately 12 000 lead glass blocks providing an azimuthal coverage up to polar angles of 
|cos#| < 0.98. To improve the shower energy and spatial resolution, a preshower detector 
(presampler) is mounted between the solenoid and the electromagnetic calorimeter. 

For background studies, about 8 million hadronic decays of the Z have been simulated using 
the JETSET 7.4 Monte Carlo || with parameters tuned to the OPAL data [IU|. In addition, 



samples with simulated D*' decays have been generated since JETSET does not include D*' 
production by default. For the production of the D*' Monte Carlo samples with JETSET, 
radially excited D meson states have been implemented into JETSET with the mass set to 
2.629 GeV/c 2 , according to the most recent calculation ||, and the width set to zero. The 
possible effects of a non-zero width are treated later as systematic uncertainty. DELPHI only 
measured the ratio of the D* /+ and Dj' 10 production rates}] 0, and therefore this measurement 

2 The OPAL coordinate system is denned as a right-handed Cartesian coordinate system, with the cc-axis 
pointing in the plane of the LEP collider towards the centre of the ring and the z-axis along the electron beam 
direction. 

3 The symbol D^*' represents the two narrow orbital resonances Di and D^. 
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has been combined with the OPAL measurement of Dj' )0 production |TTJ to calculate the 



j 

expected absolute D* /+ production rate. D*' production is simulated in fragmentation and in 
B meson decays. In the latter case, D*' mesons were produced in the same b decay channels as 
D2, in proportion chosen to reproduce the DELPHI measurement ||. 

Four samples of D*' Monte Carlo events have been generated: production in Z — > cc and 
Z — ► bb events and the decay channels D*'+ — ► D* + 7r + 7r~ and D*'+ — > D^°7r + . The signal Monte 
Carlo sample contains 600,000 events in total, with at least one charged D*' meson in each event. 
In all Monte Carlo samples, heavy quark fragmentation has been implemented using the the 
model of Peterson et al. |T2] with fragmentation parameters determined from LEP data [IB 



The samples have been passed through a detailed simulation of the OPAL detector [14] and 
then analyzed in the same manner as the data. 



3 Event and Track Selection 



Hadronic Z decays are selected based on the number of reconstructed tracks and the energy 
deposited in the electromagnetic calorimeter, as described in Ref. | 15|| . The analysis uses an 
initial sample of 3.73 million hadronic decays of the Z collected with the OPAL detector between 
1991 and 1995. 



Tracks are used in the reconstruction if they pass loose quality cuts requiring \do\ < 0.5 cm, 
\z \ < 20 cm, p > 0.5 GeV/c, and p xy > 0.150 GeV/c. Here |<i | is the distance of closest 
approach of the track to the primary vertex, measured in the plane perpendicular to the beam, 
zq the distance along the beam at this point, p the momentum, and p xy the momentum in the 
plane perpendicular to the beam. Only tracks with more than 40 hits out of a maximum of 
159 in the main tracking chamber are used. The primary vertex in a collision is reconstructed 
from the tracks in the event and constrained by the known average position and spread of the 
e + e~ interaction point. To improve the mass resolution, a track is required to have a precise 
z measurement at the exit point of the tracking chamber, either from an associated hit in the 
^-chamber surrounding the central drift chamber, or from the presampler. In cases where the 
particle exited the tracking chamber in the endcap region, i.e. | cos^l > 0.73, the exit point is 
determined precisely from the last sense wire used for the track measurement. 



4 Reconstruction of D* /+ Candidates 



The reconstruction of D* /+ — > D* + 7r + 7r decays involves D* + reconstruction and, subsequently, 
the combination of accepted D* + candidates with additional pion candidates. The D* + recon- 



struction follows closely a recent OPAL publication [J6|, with the exception of the cut on the 
kaon particle identification probability which has been tightened for this analysis, in order to 
reduce combinatorial background. 

The following decay sequence is reconstructed exclusively: 
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U D°7T + 

U K-7T+ 



All tracks that fulfil the criteria described in Section |3| are considered for the construction of 
D* /+ candidates. Only combinations of tracks that have the correct charge assignments, and 
which pass intermediate requirements imposed by the reconstruction of the D* + , are taken into 
account. The two pions produced in the initial D* /+ decay will be referred to as "D* /+ pions", 
whereas the pion coming from the D* + decay will be referred to as "slow pion" . 

Two tracks of opposite charge are accepted as a D° candidate if their invariant mass lies 
within the range 1.79 — 1.94 GeV/ c 2 , assigning the pion mass to one particle and the kaon mass 
to the other || . D° candidates are combined with tracks of charge opposite to that of the kaon 
candidate to form D* + candidates. The scaled energy xe of the D* + candidate, i.e. the ratio of 
the energy sum of all participating tracks over the beam energy, is required to be larger than 
0.2 in order to reject combinatorial background from low-energy fragmentation tracks. The 
difference of the invariant masses of D* + and D° candidates, Am = m^** — m D o, must lie within 
the interval from 142 to 149 MeV/c 2 . 

Background in the sample is further reduced by cutting on the helicity angle 9* measured 
between the direction of the D° candidate in the laboratory frame and the direction of the kaon 
in the rest frame of the D° candidate. The kaon candidate from the D° decay is expected to be 
isotropically distributed in cos^*, while the background peaks at cos6>* = — 1 and, particularly 
at small xe, at cos#* = +1. This effect is taken into account by requiring cos8* > —0.9 for 
x E > 0.5 and | cos#*| < 0.8 for x E < 0.5. 

At low xe, where the background is most pronounced, the particle identification power of 
the OPAL detector is used to enrich the sample in true kaons from D° decays. A probability 
Wk is computed from the ionization energy loss measurement of a track in the drift chamber, 
the track momentum, and the theoretical expectation for a kaon. At least 20 dE/dx out of a 
possible 159 measurement points and Wk > 10% are required for the kaon candidate in a D* + 
candidate with xe < 0.5. 

A fraction (1.5%) of all D* + candidates share the same slow pion candidate with another 
D* + candidate in the same event. In this case only the D* + candidate with a D° invariant 
mass closest to the mean value obtained from a fit to the D° mass distribution is kept for the 
further analysis. The D* + candidates surviving the cuts described above are combined with 
all remaining pairs of oppositely charged tracks. For these combinations, a mass difference 
Am = m^*/ — mo* is calculated. 

In order to suppress background from uds events, charm and bottom tags are used. The 
mean fractional energy of D* /+ mesons in Z — > cc events is expected to be large compared to 
D* /+ from B decays, and especially compared to fake D* /+ candidates in events of all flavours 
consisting of low-energy fragmentation tracks. Since most of the energy of the D* /+ is transferred 
to the D* + , a charm enriched sample is selected by imposing a cut on the energy fraction of 
the D* + of Xe > 0.4. Additionally, in the Monte Carlo simulation, the D* /+ pions have a 
higher momentum if the D*' is produced in cc events compared to the production via the 
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decay of b hadrons. This information has been exploited by selecting only candidates for which 
the magnitude of the vector sum of the two D*' pion momenta is greater than 3.6 GeV/c 
(Figure [I]). The mass resolution in the charm enriched sample is improved by constraining all 
D* /+ candidate tracks except those from the D° decay to the primary event vertex. 

A bottom-enriched sample is selected by requiring an apparent D° decay length of at least 
0.3 mm, defined as the distance between the reconstructed D° vertex and the primary event 
vertex, measured in the plane perpendicular to the beam. D° candidates from light quark (uds) 
events and cc events are expected to be made up from tracks in the vicinity of the primary 
vertex. In contrast, D° candidates in bb events tend to have a larger decay length with respect 
to the primary vertex, due to the B meson lifetime. 

The separation significance of two vertices is defined as the distance between them in the 
plane perpendicular to the beam, divided by the uncertainty on this quantity. The separation 
significance of the D° vertex and the reconstructed vertex formed by the D*' pion tracks is 
required to be between —2.0 and +4.0. The distribution of this variable is not centred at zero, 
because the D° has a lifetime large enough to be observed. The separation significance between 
the D* /+ pion vertex and the D° vertex is shown in Figure |2|. A positive sign indicates that 
the D° vertex is farther from the primary interaction point than the D* /+ pion vertex. Positive 
values are expected for signal, because the D° occurs later in the decay chain. 



flavour 


charm enriched 


bottom enriched 


combined sample 


uds 


12% 


10% 


11% 


c 


60% 


15% 


31% 


b 


28% 


75% 


58% 



Table 2: Flavour composition of the charm enriched, bottom enriched and combined Monte Carlo 
samples. 



The Monte Carlo flavour composition of the D* /+ candidates passing all selection criteria is 
shown in Table |2|. For D* /+ mesons from primary charm and bottom quarks, the reconstruction 
efficiencies are estimated from Monte Carlo analysis to be 14.1% and 7.1%, respectively. The 
combined efficiency for D* /+ from both sources is 11.4%, assuming equal rates for production 
in charm and bottom events. This assumption is motivated by the range of production rate 
ratios observed for D* ± [16] and orbitally excited Dj^° mesons [11]. The analysis has also been 
performed with several other assumptions, and the dependence of the results on these different 
scenarios is included SIS db systematic error. 



5 Selection Results 



In order to maximize the expected sensitivity of the analysis, the charm and bottom enriched 
samples with 1765 and 3051 candidates, respectively, in the Am = uld*' — hid* region between 
2.3 GeV/ c 2 and 3.0 GeV/ c 2 are merged into a combined sample. The overlap of 324 candidates 
between charm and bottom sample is taken into account; candidates that have been tagged as 
charm and bottom are used only once. This sample contains a total of 4492 D* /+ candidates 
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in 2192 events. The results for the combined sample and the b and c enriched samples are 
discussed below. 



The Am distribution of the combined charm and bottom enriched sample is shown in 
Figure |3]a. No narrow resonance is observed anywhere in the Am region between 0.3 and 
1.0 GeV/ c 2 , although the Monte Carlo simulation with a D*' production rate fixed at the value 
published by DELPHI shows a clear signal (see Figure ^b). This is also separately true for 
the c and b enriched samples, shown in Figures f| and [|, respectively. Shape and normalization 
of the background differ in data and simulation, especially in the b enriched sample. This is 
due to D* /+ candidates in the Monte Carlo where one track is incorrectly assigned to the D* /+ 
candidate. These candidates lead to a broad enhancement in the mass distribution. The mean 
scaled energy xe of these candidates is small. Therefore, they contribute more to the b enriched 
sample than to the c enriched sample. The analysis has also been performed on a Monte Carlo 
sample without D*'(2629) production. It has been found that in this case the background shape 
and normalization agree very well with the data. This absence of a significant enhancement over 
the expected background in data provides additional confidence that the D* /+ (2629) production 
rate must be small. 



A limit on the D*' + (2629) production rate is calculated by defining a mass window in the 
Am range of 0.58 GeV/c 2 < Am < 0.66 GeV/c 2 , corresponding to a D*' mass window of 
2.59 — 2.67 GeV/ c 2 . This includes the ±2cr range of both the theoretical prediction of the D*' 
p[ and the excess observed by DELPHI . The background distribution has been fitted using 
a parametrisation of the form @. I I 



f(Am') = a (Am - mof e " 7(Am ""^ , (1) 

where m® = 0.28 GeV/c 2 is the kinematic limit for D* + 7r + 7r _ combinations, and a, (3 and 
7 are the fit parameters. The signal region has been excluded from the fit. The number of 
D* /+ candidates in the mass window is determined to be 14 ± 28 by subtracting the integrated 
interpolated background function from the observed number of candidates in the signal interval. 
The quoted statistical error is obtained by combining the Gaussian variance of the number of 
entries in the mass window (±25) and the uncertainty on the background integral obtained by 
propagating the uncertainties on the fit parameters (±14), taking into account the correlations. 



6 Systematic Checks 



The D* /+ analysis has been checked by applying a similar selection to D* + 7r~ final states, looking 
for the narrow orbital resonances D^ and Dg . The D* + selection criteria are as described earlier. 
Accepted D* + candidates are combined with pion candidates passing selection criteria similar 
to those applied to the D* /+ pions, except for the cuts on the pion-pion momentum sum and 
the pion-pion vertex separation, which are inappropriate here. Instead, the pion momentum is 
required to be larger than 2.0 GeV/c, as was done in an earlier OPAL publication on orbitally 
excited D mesons [O. The efficiency of the Dj — > D* + 7r~ reconstruction is about 7.3%. The 



results achieved on data and simulated events can be compared in Figure |6|. The method used 
for determining the number of signal events is the same as was used for the D* /+ . The lower 
mass window boundary has been set to 2.382 GeV/ c 2 , two T>\ widths below the world average 
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mass. The upper boundary has been chosen to be 2.507 GeV/ c 2 , two D^ widths above the 
world average Dt, mass. The production rates of the two narrow Dj resonances have been 



adjusted in the simulation to match the previous OPAL measurement fill . The amount of 
signal found in Monte Carlo is 281 ± 16, whereas in data the excess is 189 ± 41 in a sample of 
4711 selected Dj candidates. The analysis presented here and the previous OPAL publication 
on orbitally excited D mesons [jTlJ] that was used to adjust the Dj production rates in Monte 
Carlo use a similar dataset. Both results are therefore expected to be correlated. Nevertheless, 
it has been found that the overlap of the selected Dj'* candidate samples of old and new 
analysis is very small, and the results can thus be considered almost statistically independent. 
This leads to the conclusion that the measured rates are statistically compatible. 

Having confirmed in the case of D* + 7r~ final states that the reconstruction procedure per- 
forms as expected, the results of the search in the D* + 7r + 7r~ final state are used to place limits 
on the production rate of the D* /+ (2629). The efficiency of D* /+ reconstruction is taken from 
Monte Carlo simulation. The efficiency calculation accounts for several sources of systematic 
uncertainties. They are evaluated as follows: 



The efficiency has been calculated from a simulation with zero D* /+ width. To quantify the 
effect of a non-zero width, the analysis has been repeated with a D* /+ width of 15 MeV/ c 2 , 
the 95% C.L. upper limit obtained by DELPHI |5]]. This results in a relative decrease of 
the efficiency of 16.5%, taken as systematic uncertainty. 

The simulated detector resolutions of momentum, impact parameters and track angles 



have been varied by ±10%, as in previous publications [II], |IE , around the values that 

h4.9?^ 
-5.7%/ 



best describe the data. The largest relative deviation in the selection efficiency (ti'ly) is 



taken as systematic uncertainty. 

Equal production rates of D* /+ in c and b decays have been assumed. Other possible 
assumptions include the same ratio as observed in D* 1 * 1 production [[HJ, or in Dj pro- 



duction |il |. Equal numbers of D*' from both primary heavy flavours are simulated as 
well as the ratio of candidates from bb and cc events observed by DELPHI ||, corrected 
for the different reconstruction efficiencies of the respective analyses. The ratios thus 
obtained cover a 1^3% ra nge around the central assumption of equal production rates in 
bb and cc events. The largest observed variations in the efficiency in both directions, 
^46%, are taken as systematic uncertainty. 

A Peterson function has been used to parametrise the heavy quark fragmentation in 



the Monte Carlo samples. Other fragmentation models by Collins and Spiller |L7j and 
Kartvelishvili |18] have also been tested. The fragmentation parameters have been varied 
in the limits implied by the measurement of the mean xe distribution of D* + mesons at 
LEP [13| . The largest resulting uncertainty of ^ 3% * s taken as systematic uncertainty. 

The efficiency depends on the performance of the particle identification by dE/dx mea- 
surements. The dE/dx calibration for kaons has been checked in a previous analysis [|l]| 
under identical circumstances. An error of ±3.2% was found for the total rate of kaons 
passing the selection. Because this cut is applied only for candidates with x^ < 0.5, the 
expected contribution to the uncertainty on the efficiency is ±2.1%. 
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• The D° lifetime is currently measured with a precision of 0.7%, whereas the average b 
hadron lifetime is known to 0.9% ||. The cut on the D° apparent decay length and thus 
the b flavour enrichment is sensitive to the modelling of this quantity. The cut value 
has therefore been varied within the combined uncertainty on D° and B lifetime. The 
resulting deviation of 1q o% * s taken as systematic uncertainty. 

• The calculated production limit depends on the branching ratio Br(D* + — > D°7r + ) x 
Br(D° -> K-7T+). The world average value is (0.677 ± 0.005) x (0.0383 ± 0.0009) @. The 
relative uncertainty on the product branching ratio is ±2.5%. 

• The reconstruction efficiencies of the decay chains D* /+ — > D^°7r + , — > D* + 7r _ and 
D* /+ — > D* + 7r + 7r _ are identical within the statistical uncertainties. No systematic un- 
certainty is introduced. 

• The reconstructed number of D* /+ in Monte Carlo has been compared to the actual 
number of D* /+ in the sample. The excess over the background fit is 138 ± 30, with 175 
entries due to D* /+ with all tracks reconstructed. It is thus concluded that the fit is free 
of significant bias. 

• The width of the signal region has been varied from 90% to 130% of the nominal value 
to check whether the resulting counting rate variations are consistent with statistical 
effects. All results are clearly within the statistical error. Thus no additional systematic 
uncertainty is assigned. 

• The bin width of the Am histograms in which the fit is performed has been varied from 
75% to 125% of the value used for the analysis. No significant impact on the rate is 
observed. 

The contributions to the systematic uncertainty are summarized in Table |[ 



7 Calculation of D* /=t (2629) Production Limits 

The limit is calculated assuming that the candidate sample in the signal region is composed 
of a large background and a small signal sample. Each contribution has a Poisson probability 
density function although the background can be approximated by a Gaussian distribution 
ebck, o) since the background expectation is large. The width of the Gaussian, a = 28, has 
two contributions, one accounting for the statistical uncertainty in the signal window, the other 
one taking into account the fit error. n Q = 610 candidates are observed in the signal window. 
The probability to count no candidates or less is given by 

no n 

P(n ;e sig ) = J2 J2 V(n sig ;e sig ) g(n-n sig ;e hck ,a) (2) 

n=0 n s i g =0 

where e S i g is the unknown expectation value of the signal and V its Poisson distribution. The 
summation is performed over all possible combinations of n, the total number of candidates 
in the mass window, and n S j g , the amount of signal within these n candidates, where n < Uq. 
Assuming no narrow D* /+ are produced, i.e. e S i g = 0, a probability of 0.70 is obtained to observe 
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error source 


relative contribution 




c enriched 


b enriched 


combined 


relative errors on Monte Carlo efficiency 


D*'+ width 


+0.0% 
-17.9% 


+0.0% 
-10.6% 


+0.0% 
-16.5% 


detector resolution 


+5.2% 
-6.1% 


+4.5% 
-5.7% 


+4.9% 
-5.7% 


relative production rates in b and c 


none 


none 


+3.8% 
-4.6% 


fragmentation modelling 


+3.6% 

-3.7% 


+1.1% 

-0.3% 


+1.7% 

-2.3% 


kaon dE/dx 


+1.3% 

-1.3% 


+2.6% 
-2.6% 


+2.1% 
-2.1% 


B and D° lifetimes 


none 


+0.5% 
-0.0% 


+0.2% 
-0.0% 


Monte Carlo statistics 


+0.1% 
-0.1% 


+0.08% 
-0.08% 


+0.07% 
-0.07% 


total 


+6.5% 
-19.3% 


+5.3% 
-12.3% 


+6.8% 
-18.4% 


relative errors on external branching ratios 


branching ratio D* + — > D°7r + , D° — > K rr + 


±2.5% 


±2.5% 


±2.5% 


error on r bE /r had and r c c/r had 


±2.9% 


±0.3% 


none 


total 


+7.5% 
-19.7% 


+5.9% 
-12.6% 


+7.2% 
-18.6% 



Table 3: Overview of the systematic error sources contributing to the total uncertainty. 



610 candidates or less. Given the prior knowledge of the maximum amount of background, the 
95% C.L. limit is obtained at P(n ; e sig ) = 0.70 x 0.050. At large n 10) the systematic 
uncertainties affecting the efficiency start playing a role in calculating the limit. They are 
considered by substituting the Poisson distribution of the signal with a Gaussian distribution 
at n > 10 and adding in quadrature the asymmetric systematic uncertainty of the efficiency 
to the width of the Gaussian given by y/n. The obtained limit on the production of narrow 
D*'+.(2629) j n hadronic Z decays is: 

/(Z -> 0*^(2629)) x Br(D* /+ -> D* + 7r+7r-) < 3.1 x 10~ 3 (95% C.L.). (3) 

This corresponds to a 95% C.L. upper limit of 66 on the number of reconstructed D* /+ in 
the signal region. The approach used here is similar to a Bayesian approach with a flat prior 
distribution for positive D* /+ production rates and zero elsewhere, but superior in the fact that 
negative expectation values are excluded in principle by using the proper Poisson distribution. 

Limits on the D* /+ (2629) production rate in charm and bottom events are also calcu- 
lated. In the charm and bottom enriched samples, the numbers of events in the signal re- 
gion, relative to the expectations derived from the fitted background functions, are 5 ± 18 
and 29 ± 23, respectively. In using these results to calculate limits, a further systematic 
uncertainty arises due to the experimental uncertainties on r c c/Phad = 0.1671 ± 0.0048 and 
Lkb/Lhad = 0.21644 ± 0.00075 0. Under the conservative assumption that D* /+ are only pro- 
duced in Z — > cc events, i.e. that the measured excess can be fully assigned to D*' + candidates 
in cc events, a production limit of 

/(c -> D* /+ (2629)) x Br(D*'+ -> D* + 7r+7r) < 0.9 x 10~ 2 (95% C.L.) (4) 
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is obtained from the charm enriched sample, while under the opposite assumption that all D* /+ 
are only produced in the decay of b-hadrons, a limit of 

/(b -> D*'+(2629)) x Br(D*'+ -> D* + 7r + 7r-) < 2.4 x 1(T 2 (95% C.L.) (5) 

is computed using the bottom enriched sample. The systematic uncertainties assumed for the 
calculations in the separate samples are given in Table |3|. 



8 Discussion 



The result of the OPAL D* /+ search does not agree with the published DELPHI results ||. 
Both analyses apply standard selection criteria to obtain a high-purity D* + sample. In contrast 
to the OPAL analysis, DELPHI reconstructs D° candidates in two decay channels, D° — > K~7r + 
and D° — > K~7r + 7r + 7r~. The loss in efficiency associated with using only the K~ir + channel is 
in part compensated by the softer kinematic cuts and a slightly larger hadronic event sample. 
Combinatorial background is reduced by applying a best candidate selection to D* + candidates 
which share tracks. These candidates can produce correlated entries in the Am distribution 
if they are combined with the same D* /+ pions. The expected sensitivity of the analysis to 
the existence of a narrow resonance in the mass region including the theoretical predictions for 
the D*' as well as the mass of the published DELPHI result, produced at a rate published by 
DELPHI, is demonstrated in Figure |^. The excess observed in the Monte Carlo sample with 
comparable size to the data sample has a significance of 4.7a (statistical error only), or 3.5<r 
including all systematic uncertainties. 

The OPAL and DELPHI selection criteria used to select charm and bottom enriched sam- 
ples, and for the pions from the D* /+ decay, have been optimized in different ways. For example, 
the kinematic cuts for the pion candidates from the D* /+ are less restrictive, 0.5 GeV/ c in this 
analysis while DELPHI requires 1.0 GeV/c. 

The sensitivity of the OPAL analysis has been checked with the Dj reconstruction de- 
scribed above, where agreement between data and simulation is observed. The selection criteria 
that are specific to either the D* /+ or analysis are checked by performing the analysis with- 
out these cuts. This did not change the result of the D* /+ search. 

Overall, the sensitivities of the two analyses are found to be similar, although the background 
in this analysis is higher. A comparison between the analyses can be made by calculating the 
D* /+ versus Dj*' )0 rate for which DELPHI published a number M: 

R „ (AW)Br(D*'+ -> D*+7r+7r-) 

(iV D o)Br(D? -> D*+tt-) + (iV D *o)Br(Df -> D*+tt-) 1 ' 

= 0.49 ± 0.18(stat.) ± O.lO(syst) , 

where (Nd*/+) denotes the expected number of D* /+ in a sample with (N D o) and (iV D *o) being 
the corresponding quantities for the D° and D^ . Based on the numbers presented in Sections 
|] and H, 

R = 0.05 ± 0.10(stat.) ± 0.002(syst.) (7) 
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is calculated. The statistical error of the fit and the systematic uncertainties resulting from 
the D*' width and relative production rates in b and c are included. In the ratio calculation, 
all other systematic uncertainties largely cancel and have therefore been neglected. Using a 
statistical approach analogous to the one described earlier, a limit of 

R < 0.22 (95% C.L.) (8) 

is computed. 

In summary, the evidence of D* /+ production in hadronic decays of the Z published by the 
DELPHI collaboration |J is not confirmed with OPAL data. 
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Figure 1: Magnitude of the vector sum of the D*' pion momenta, p n7T , for data (points with error 
bars), Monte Carlo without D* /+ (open histogram, scaled to the same number of entries as data) and 
true D* /+ (shaded, scaled up by a factor of 10). The arrow indicates the selected region. 
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Figure 2: Separation significance a between the D° vertex and the reconstructed vertex of the D* /+ 
decay pion tracks for data (points with error bars), Monte Carlo (open histogram, scaled to the same 
number of entries as data) and true D* /+ (shaded, scaled up by a factor of 10) for candidates passing 
all selection criteria except for the cut on the quantity shown here. The arrows indicate the selected 
region. 
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Figure 3: Distribution of the mass difference Am = ran" — m^* o/D* /+ candidates in the combined 
charm and bottom sample, for (a) data and (b) Monte Carlo. The calculation of Am as well as 
the background parametrisation superimposed on the histogram are described in the text. The Monte 
Carlo histogram is scaled to the number of hadronic events in data. Furthermore, the D* + Monte 
Carlo production rate has been adjusted to an OPAL measurement of this quantity fijly, and the D* /+ 
production rate has been fixed at the value published by DELPHI and the OPAL T)j° measurement 



ftlJ\] . The Monte Carlo plot presented here is created from a subsample of the available events that has 
roughly the same size as the data sample. The arrows indicate the mass window defined as the signal 
region. The line represents the result of the background fit, where the line is dashed in regions that 
have been excluded from the fit. The X 1 /dof of the background fit is 1.23 for data, and 0.99 for the 
Monte Carlo simulation. 
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Figure 4: Distribution of the mass difference Am = m^./ — mc* of D* /+ candidates in the charm 
enriched sample, for (a) data and (b) Monte Carlo. The X 2 /dof of the background fit is 0.80 for data, 
and 0. 70 for the Monte Carlo simulation. 
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Figure 5: Distribution of the mass difference Am = m^*/ — mo* of D* /+ candidates in the bottom 
enriched sample, for (a) data and (b) Monte Carlo. The x 2 /dof of the background fit is 1.20 for data, 
and 1.07 for the Monte Carlo simulation. The difference in shape and normalization of the background 
is explained by a large number of partially reconstructed D* /+ candidates which, due to their low mean 
scaled energy, mainly contribute to the b enriched sample. The analysis has also been performed on a 
Monte Carlo sample without D*'(2629) production. It has been found that in this case the background 
shape and normalization agree very well with the data. 
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Figure 6: Distribution of the mass Am„ = m D *+ 7r - — m D *+, /or (aj data and (b) Monte Carlo. 
TTie production rates of D j mesons in the Monte Carlo have been adjusted to the values measured 
in data J71|/ . The same parametrisation as in the D* /+ analysis is used to describe the background. 
The x 2 /dof of the fit is 1.43 for data, and 1.22 for the Monte Carlo simulation. The Monte Carlo 
histogram is scaled to the number of hadronic events in data. 
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